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Summary 

An analysis has been performed to compare de> 
coupled and linear quadratic regulator (LQR) pro- 
cedures for the control of a large, flexible space an- 
tenna. Control objectives involved (1) commanding 
changes in the rigid-body modes, (2) nulling initial 
disturbances in the rigid-body modes, or (3) nulling 
initial disturbances in the first three flexible modes. 
Control was achieved with two three-axis control- 
moment gyros located on the antenna column. The 
results are presented to illustrate various effects on 
control requirements for the two procedures. These 
effects include errors in the initial estimates of state 
variables, variations in the type, number, and loca- 
tion of sensors, and deletions of state-variable esti- 
mates for certain flexible modes after control acti- 
vation. The advantages of incorporating a time lag 
in the control feedback are also illustrated. In addi- 
tion, the effects of inoperative-control situations were 
analyzed with regard to control requirements and re- 
sultant modal responses. Comparisons are included 
which show the effects of perfect state feedback with 
no residual modes (ideal case). Time-history re- 
sponses are presented to illustrate the various effects 
on the control procedures. 

Introduction 

The advent of the NASA Space Transportation 
System (STS) has provided the capability of estab- 
lishing large space structures in Earth orbit for sci- 
entific and operational missions. Such structures in- 
clude large, flexible space antennas and the proposed 
space station, which may be, in part, flexible. These 
structures will require control for attitude orienta- 
tion as well as procedures for controlling the flexible 
modes. These control problems are considered in this 
paper for a 122-m (400-ft) diameter hoop-column an- 
tenna (ref. 1). A finite-element mathematical model 
of the antenna is used, which provides an excellent 
tool for simulated control analysis of realistic large, 
flexible space structures. Previous control studies of 
this antenna model include references 2 and 3. The 
antenna is depicted in sketch A and is described in 
reference 1. 

The present paper is an extension of the study 
reported in reference 2. This study represented the 
idealized condition; that is, perfect knowledge of the 
state vector was assumed for feedback and no residual 
modes were considered in the mathematical model. 
Most of the results pertained to a decoupled-control 
procedure with a brief comparison shown for a lin- 
ear quadratic regulator (LQR) procedure. The re- 
sults were presented essentially as a parametric anal- 
ysis of rigid-body and flexible-mode control require- 
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Solar panels 


Sketch A 

ments over a large range of closed-loop frequencies 
and damping ratios. Various arrangements and lo- 
cations for two types of control actuators (control- 
moment gyros (CMG’s) and reaction-control jets) 
were investigated. Effects of varying the number of 
controlled flexible modes were also determined. 

In the present paper, an observer is used for es- 
timating the state variables. In addition, residual 
modes are incorporated into the analysis and their 
effects are investigated. These modes can be as- 
sumed to be either unknown (not included in the 
antenna model) or known but not included in the 
control law. In reference 2, the control-actuator ar- 
rangement which appeared to be the most practical 
employed one three-axis CMG at the top of the col- 
umn and one at the bottom. With one exception, this 
arrangement is used for the present analysis. The val- 
ues selected for the closed-loop dynamics are based 
on the results of reference 1. Only the first three 
flexible modes are included in the control law, be- 
cause the data in reference 2 showed large increases 
in the control requirements with additional modes 
included. It was also shown in reference 2 that the 
decoupled- and LQR-control procedures gave com- 
parable results; hence, both procedures are analyzed 
and compared in the present analysis. 

Effects on the control requirements are analyzed 
herein for (1) commanding changes in the rigid-body 
modes, (2) nulling initial disturbances in the rigid- 
body modes, or (3) nulling initial disturbances in the 
first three flexible modes. The study includes the ef- 
fects of (1) inoperative-control situations, (2) values 
selected for closed-loop dynamics, (3) errors in ini- 
tial knowledge of state variables, (4) time lag in con- 
trol feedback, (5) number and location of sensors, 
and (6) deletions of estimates in some of the flexible 
modes. Results are shown in tabular form, in para- 
metric plots, and as sample time histories of modal- 



amplitude and control responses. The control sys- 
tem and observer are assumed to be perfect, in that 
no actuator or sensor dynamics are included in the 
analysis. 

Symbols and Abbreviations 

A system matrix (eq. (4)) 

An modal amplitude (eq. (2)), where n 

is the mode number (1,2,..., 6) 

B control influence matrix (eq. (4)) 

C observation matrix (eq. (9)) 

CMG control-moment gyro 

d displacement vector (eq. (3)) 

F decoupled feedback gain matrix 

(eq. (5)) 

G decoupled feedforward gain matrix 

(eq. (5)) 

I moment-of- inertia matrix (eq. (1)) 

Ix^Iy.Izy Ixz ceuter-of- gravity moments and 
product of inertia 

J objective function (eq. (7)) 

K LQR feedback gain matrix (eq. (6)) 

K' estimator gain matrix (eq. (13)) 

LQR linear quadratic regulator 

I distance (see fig. 1) 

Mij angular momentum, where i is 

the control direction and j is the 
actuator location 

run modal mass 

P control-vector weighting matrix 

(eq. (7)) 

Q state-vector weighting matrix 

(eq. (7)) 

r radius of hoop 

s vector of system input noise 

(eq. (10)) 

T torque vector 

Ti^j components of torque, where i is 

the control direction and j is the 
actuator location 

Tmax maximum value for actuator requir- 

ing largest torque 


Tn natural period 

r^otal sum of maximum torque values for 

the six actuators 

t time 

u observation noise (eq. (12)) 

V input command vector (eq. (5)) 

w weight 

z coordinates of antenna center of 

gravity 

X state vector 

y observation vector (eq. (8)) 

z sensor output vector without noise 

s estimator error, y — C'x 

desired damping ratio 
natural damping ratio 

9^(j),'il) rotation angle about a;-, y-, and 0 - 

axis, respectively 

T time constant 

j mode-shape matrix, where n is the 

mode number, i is the direction, 
and j is the location 

f(<n) 

j ^ mode-slope matrix, where n is the 

mode number, i is the direction, 
and j is the location 

(jJd desired frequency 

ujn natural frequency 

Subscripts: 

F feedback 

G feedforward 

control direction and actuator 
location, respectively 

A matrix with a prime relates to the estimator 
equations. Dots over symbols indicate derivatives 
with respect to time. A circumflex ( " ) over a symbol 
indicates an estimate of the state variable. All 
matrices are given in units of inches, pounds, and 
seconds. 

Mathematical Model of Antenna 

A large space structure such as the hoop-column 
antenna has, in theory, an infinite number of flexi- 
ble (vibration) modes. To facilitate analytical treat- 
ment of the control problem, a finite-order, linearized 
model was formulated. For the analysis of this re- 
port, the structural model was selected to contain 
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six of the lowest flexible modes of the 122-m-diameter 
hoop-column antenna described in reference 1 . The 
three rigid-body rotation modes are included in the 
analysis. 

The equations of motion used to represent the 
rigid-body and vibration modes of the antenna are 
given below. Rigid-body rotations (for small angles) 
about the antenna center of gravity for two three-axis 
CMG actuators on the column are represented by 


where the state vector 

x= [ 6 , (/>, Ai, A2, . . . , Ae, A2, Ae]^ 

Ai, A2, and A3 are the controlled flexible modes, and 
A4, A5, and Ag are considered as residual modes. 
The B matrices used in the present analysis are 
given in table II. The composition of the A matrix is 
deflned in reference 2. 


where T — ^2/ A’ ^^, 1 ’ ^ The 

CMG’s were used on the column such that control 
torques were about the x-, y-, and ^-axes, as defined 
in figure 1. The CMG’s are numbered to be consis- 
tent with reference 2. 

Variations in modal amplitudes of the flexible 
modes are represented by 


rrin 


(n-l,2,...,6) (2) 

where rrin is the modal mass and is the mode-slope 
matrix. It should be noted that values in equa- 
tion (2) are modal-amplitude displacement variables 
and do not represent actual physical displacements. 
The physical displacement at some point on the an- 
tenna includes linear combinations of the modal am- 
plitudes and the mode shapes and is given by the 
transformation 

dj = $An ( 3 ) 

where ^ is the mode-shape matrix. 

The mode-slope data used in the analysis wens 
taken from unpublished results of a NASTRAN® 
model of the 122-m-diameter hoop-column antenna 
system and are provided in reference 2. Table I 
shows the weight and the inertias of the antenna. 
The modal masses, the natural frequencies, and the 
natural damping ratios of the six flexible modes 
considered in the analysis are also given in table I. As 
a matter of interest, the unforced transient behavior 
of these modes is presented in figure 2. 

Decoupled and LQR Control 

The second-order equations in the analysis 
(eqs. (1) and (2)) can be reduced to first-order equa- 
tions (state-vector form) and written as 


X = Ax + BT 


( 4 ) 


Decoupled Control 

Since decoupling theory is analyzed extensively 
in the literature, the discussion of its application in 
this paper is limited to some general remarks. The 
decoupling control law is taken as 

T = Px -}- Gv (5) 

where x is the estimated state vector, v is the in- 
put command vector, and F and G are feedback and 
feedforward gain matrices, respectively. (As noted, 
estimates of the modal variables are used because 
these variables cannot be measured directly.) De- 
coupling theory provides a method for determining 
the F and G matrices such that independent control 
is maintained for each of the decoupled (controlled) 
state variables. Different values can be selected for 
the closed-loop dynamics (w^ and f^/) without affect- 
ing the independent control capability. Examples are 
given in reference 2 which illustrate the manner in 
which the theory was applied to the system equa- 
tions in determining the F and G matrices. The 
values for the closed-loop dynamics used in the cur- 
rent analysis are given in table III. The “standard” 
case corresponds to Q values which give approxi- 
mately the same value (10 sec) for time to damp 
to 1 percent for each of the three -flexible modes. 
The corresponding F matrix for the standard case 
is shown in table IV for two control arrangements. 
The G matrix is given in table V. An additional 
F matrix, used for LQR comparison, is given in ta- 
ble VI. The F and G matrices were calculated for a 
reduced-order model (no residual modes); hence, x = 

[^, Ai, A2, A3, d, 0 , Ai, A2? A3]^. When 
the residual modes are included in the estimator, 
zero columns must be added to the F matrices in 
order to make the matrices of equation (5) compati- 
ble. This has been done in table IV by adding three 
zero columns after both column 6 and column 12 . In 
this case, x becomes an 18 x 1 vector. The decoupled 
and LQR gain matrices, as well as all other matrices 
in the paper, are given in units of inches, pounds, 
and seconds. 
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LQR Control 

The asymptotic LQR approach of reference 4 was 
also used to calculate the state feedback control 
gains. (A computer program is given in ref. 5 for 
performing these computations.) The control law for 
the LQR results is given by 

T = Kx (6) 

The control gain matrix is computed subject to 
the constraint of equation (4) such that the following 
function is minimized: 

rt 

J = (x^Qx + T^PT) dt (7) 

where Q and P are positive definite symmetric 
weighting matrices for the states and controls and 
can be varied in order to achieve desired closed-loop 
dynamics. For the present analysis, the identity ma- 
trix was used for P; the state- vector weighting ma- 

Incorporation of Observer 


trix Q was adjusted to obtain three different sets 
(cases) of closed-loop dynamics. (See table III.) 

The K gains were used in the control law (eq. (6)) 
for the subsequent LQR control analysis. In addition, 
in order to compare the decoupled- and LQR-control 
procedures, the resultant LQR closed-loop dynamics 
(cu^ and <;^ were used as inputs to the decoupling 
program to calculate the decoupled gains. Most of 
the comparison results pertain to case 2 (table III); 
hence, the gain matrices for the two control proce- 
dures (both having identical closed-loop dynamics) 
are presented in table VI for this case. The weight- 
ings Q used for the LQR procedure are given in ref- 
erence 2. The B matrix is the same as that shown at 
the top of table II. The K gain matrices for cases 1 
and 3 are given in reference 2. (With the observer 
incorporated into the system, the 6 x 12 matrix was 
expanded by adding three zero columns after both 
column 6 and column 12 to account for the three 
residual modes.) 


An observer was incorporated into the control procedures in order to calculate estimates of the state 
variables as required in equations (5) and (6). An attitude-measurement device (e.g., a star tracker) and a 
rate-measurement device (e.g., a rate gyro) were assumed to be present at certain locations on the column. 
The sensor output is given by 

y = z + u (8) 

where z — Cx. As an example, the observation matrix C is 
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(9) 


where 13 x 3 is the identity matrix and z = [0, (j>^ ip, Ai, A 2 , . . . , Ag, 0, Ai, A 2 , . . . , Ae]^. The super- 
scripts on the mode slopes refer to the flexible-mode number. The subscripts on the mode slopes refer to the 
direction and location. For example, the subscript a:, 1 refers to the mode slope in the ^-direction at location 1 
on the column (fig. 1). The observation matrix shown in equation (9) senses attitude rate in the a:-, y-^ and 
2 f-directions at location 1 and senses attitude in the x-, j/-, and ^:-directions at location 4. Mode-slope values 
for various locations on the column are given in reference 2. The values corresponding to equation (9) are 
presented in table VII. 
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The complete set of equations required for the control analysis becomes 


X = Ax + BT + s (10) 

T = Kx (for LQR control) (Ha) 

T = Fx + Gv (for decoupled control) (Hb) 

y = Cx + u (12) 

^ = A'x + B'T + K'(y - C'x) 

= (A' - K'C' + B'K)x + K'Cx (for LQR control) (13a) 

= (A' — + B'F)x + K'Cx + B'Gv (for decoupled control) (13b) 


The estimator equations (eqs. (13)) calculate the estimates of the states jc from the measurement equation y 
and the estimator gains K', which are determined with a full-order, steady-state Kalman-Bucy filter. If all 
modes are estimated, C' = C. However, if, for example, the residual modes are not estimated, C' becomes a 
6 X 12 matrix. The estimator gains were calculated with the technique given in reference 
The set of composite equations used for the LQR-control analysis is 


•^18x18 I Bi8x6 K6x12 


LK^12x 6 CexlS I A'i2x12 — K'i2x6 C'6x12 + B'i2x6 1^6x12 
For the decoupled-control analysis, K is replaced by F in equation (14) and the term 

Bi8x6 Gqx6 

V 

,B^12x6 (*6x6. 

is added to the right side of the equation. The subscripts indicate the dimensions of the various elements for 
the case in which the three residual modes are not estimated. With precalculated values for K' and K (or K', 
F, and G), the composite set of equations was integrated to obtain time-history responses for specified initial 
conditions in the various modes. In determining the estimator gain matrix, the weighting matrices S and U for 
the states and the observations, which are analogous to the covariances of s and u, were varied until acceptable 
performance (satisfactory estimator eigenvalues and realistic gain values) was obtained. The resulting matrices 
were 

S = diag(.l, .1, .1, 10^, 10®, 10®, 10“, 10^, 10®, .1, .1, . . . , .1) 

U = diag(l, 1, 1, 1, 1, 10) 

for the case in which the three residual modes are estimated, and 

S = diag(.l, .1, .1, 10®, 10®, 10®, .1, .1, . . . , .1) 

U = diag(l, 1, 1, 1, 1, 1000) 

for the case in which the three residual modes are not estimated. 

A block diagram for decoupled control is shown in sketch B. (A block diagram for LQR control is the same, 
with the exclusion of the G matrix.) 
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Sketch B 


Results and Discussion 


Results of the study are presented to illustrate 
various effects on the control requirements. The ef- 
fects are analyzed for (1) commanding changes in 
the rigid-body modes, (2) nulling initial disturbances 
in the rigid-body modes, or (3) nulling initial dis- 
turbances in the first three flexible modes. The 
cases for (1) and (2) are referred to as rigid-body 
command and rigid-body control, respectively, even 
though control is being applied to maintain the first 
three flexible modes at zero. Likewise, the cases for 
(3) are referred to as flexible-mode control. Com- 
parisons are made between the decoupled- and LQR- 
control procedures, with the results shown in tabular 
form, parametric plots, and sample time histories of 
modal-amplitude and control responses. All time his- 
tories are for the observer in the control-system loop. 

Momentum and the maximum value for torque 
input are used as measures of control requirements 
when making various comparisons. The momentum 
time histories represent the area under the torque- 
response curves. The momentum values in the tables 
and parametric plots are the values at the end of 
each computer run, when the torque responses have 
been essentially zeroed out. The total torque is 
the sum of the maximum absolute values for the 
six actuators. Maximum torque is the value for 
that actuator requiring the highest absolute value of 
torque. 


General Considerations 

For simulation, two three-axis CMG’s were used 
as actuators to control initial disturbances of 0.01 rad 
in the rigid-body modes or of 1 in. in the three flexi- 
ble modes. Since the system equations are linear, the 
results can be scaled to any constant positive multi- 
ple of these initial conditions. Initial disturbances 
were also included in the three residual modes, even 
though it was determined that the results were not 
materially affected by their inclusion. Except where 
noted, the CMG’s were placed at the top and bot- 
tom of the antenna column at locations 1 and 2 as 
shown in figure 1. It was assumed that three-axis 
attitude sensors and three-axis attitude-rate sensors 
were available at various locations on the column for 
determining estimates of the state variables. Except 
where noted, one three-axis attitude sensor at lo- 
cation 1 was used for rigid-body control, and one 
three-axis attitude-rate sensor at location 1 and one 
three-axis attitude sensor at location 4 were used for 
flexible-mode control. 

In general, it was assumed that the control sys- 
tem was not activated until the initial estimates were 
known to 90 percent of the true values for the rigid- 
body modes and for the first three flexible modes. 
The initial estimates are actually the initial con- 
ditions in the estimator equations. Except where 
noted, estimates of the residual modes were included 
in the control procedure. No initial knowledge was 
assumed for the residual modes, inasmuch as these 
modes were not included in the control law. It was 
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found that the initial estimate in these modes was of 
little significance to the control requirements. 

Various sets of values were used for the closed- 
loop dynamics and The standard case 

refers to those values listed in table III, in which 
the damping ratios for the three flexible modes were 
selected such that the time required to damp the 
three flexible modes to 1 percent was approximately 
10 sec. The LQR closed-loop values are included 
in table III. The values for case 2 are close to the 
standard values. As stated previously, these values 
were used in the decoupling procedure for most of 
the comparisons with the LQR procedure. 

EflFect of Observer 

In reference 2, the effect of an observer was not 
studied; that is, perfect knowledge of the state vari- 
ables was assumed for feedback. Also, the effect of 
residual modes was not included. The present pa- 
per is an extension to the work of this reference and 
includes both effects. Results with and without an 
observer for three cases of closed-loop dynamics are 
presented in table VIII. The three cases correspond 
to the closed-loop values in table III, in which it is 
shown that the three sets of values are essentially 
the same except for the flexible-mode damping ra- 
tios, which represent high, medium, and low damping 
cases. The procedure used in arriving at these values 
is described in the section entitled “LQR Control.” 
It should be noted that the same closed-loop eigen- 
values for the decoupled and LQR controls were used 
to obtain the results given in table VIII. Note that 
the closed-loop flexible- mode frequencies in table III 
are close to the natural frequencies (table I). 

The data in table VIII show that the control re- 
quirements are appreciably higher when the observer 
is used for decoupled or LQR rigid-body control than 
when the observer is not used. For flexible-mode con- 
trol, the requirements are generally lower with the 
observer than they are without. Typical time his- 
tories of control without an observer are shown in 
reference 2, and time histories of control with an ob- 
server are shown and discussed subsequently in this 
report. 

In comparing the three cases in table VIII, the 
requirements for rigid-body control are not affected 
by flexible-mode closed-loop damping. The require- 
ments for flexible-mode control, however, decrease 
as the damping is lowered. Also, the momentum re- 
quirements are lower for LQR control than for de- 
coupled control, with a pronounced reduction shown 
for low damping (case 3). 

For case 2, data are also included for flexible- 
mode control with one three-axis attitude sensor at 


location 1. The results show that one sensor is 
sufficient, especially for LQR control. 

Transient Effects 

For ideal (no observer) decoupled control, there 
are no transients in the undisturbed modes. With 
the incorporation of an observer, transients are in- 
troduced, and they are compared in this section with 
LQR transients, which are present in either case. 

As far as control requirements are concerned, the 
LQR procedure appears to be the better, especially 
for flexible-mode control. However, LQR control has 
more of an effect on the modes that are not initially 
disturbed, that is, the modes that are intended to 
remain at zero. These transient effects, with the ob- 
server in the system, are presented in table IX and 
in figures 3 and 4. The data in the table show that 
the flexible-mode transients for cases 1 and 2 are 2 
to 6 times higher for LQR control than for decou- 
pled control. For low damping (case 3), these tran- 
sients are more than an order of magnitude higher. 
The flexible-mode transients increase with reduced 
damping for LQR control and decrease with reduced 
damping for decoupled control. In all three cases, 
the rigid-body transients are about an order of mag- 
nitude higher for LQR control than for decoupled 
control. For decoupled control, these transients are 
not affected by damping; however, for LQR con- 
trol, they decrease substantially as the damping is 
reduced. The transient effects may be important 
to the overall performance of the control system. 
Examples of the transients are shown in figures 3 
and 4. 

Figure 3 illustrates the transient response in the 
third flexible mode for rigid-body control. Results 
were similar for the first and second modes. Results 
are shown in which a time lag has been introduced 
into the control feedback by multiplying the feedback 
gain matrix by 1— An important characteristic 
of decoupled control is that columns of the feedback 
(or feedforward) gain matrix can be altered without 
affecting the decoupling process. The desired dynam- 
ics will be changed, but the system will still be de- 
coupled. For example, multiplying a column corre- 
sponding to a particular rate feedback by a constant 
changes the desired damping ratio by that constant. 
Columns of the gain matrix can also be multiplied by 
time- varying quantities without affecting the decou- 
pling. The overall effect is that, as time increases, the 
damping ratios increase from their natural values to 
the closed-loop values, the flexible-mode frequencies 
change from their natural frequencies to the closed- 
loop values, and the rigid-body frequencies increase 
from zero to the closed-loop values. Introducing lag 
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into the LQR-control feedback matrix gives a similar 
effect. (The stability of the system with lag added 
cannot be analytically guaranteed; however, the sys- 
tem was always stable for the magnitudes of lag used 
in the present analysis.) Additional data concern- 
ing the effect of lag on the control requirements are 
shown subsequently. 

Examination of figure 3 shows that the small 
transient in A 3 for decoupled control is virtually 
eliminated with the inclusion of lag in the control 
feedback; the transient for LQR control is reduced 
by about 50 percent. Similar effects are obtained in 
the rigid-body transients for flexible-mode control, 
shown in figure 4. In this case, lag reduces the 
transient for LQR control by about 75 percent. 

Effects on Decoupled Control Requirements 

Table X contains various effects of sensor com- 
binations on decoupled-control requirements. These 
effects include type and number of sensors as well 
as their location. Also, the effects of using differ- 
ent conditions for estimating the flexible modes are 
presented. 

The standard values for the closed-loop dynamics 
were used to generate the data in table X. For these 
data, initial estimates were assumed to have been 
determined prior to turning on the control system, 
as described previously. The different conditions 
for estimating the flexible modes, as listed in the 
columns of table X, represent conditions after the 
control system has been activated. 

The upper half of table X corresponds to rigid- 
body control and the lower half to flexible-mode 
control. Values selected for the initial disturbances 
in the residual modes have negligible effect on the 
overall results. (The values shown in the table were 
determined from data given in table IV of ref. 2 .) A 
scan of table X shows that the type, number, and 
location of sensors have, for the stable cases, only a 
small effect on rigid-body control requirements, with 
differences between the cases being about 10 percent. 
The effect on flexible-mode control requirements is 
larger, with differences being about 20 percent. The 
use of two three-axis CMG’s produces the lowest 
control requirements. 

In table X, the data in the first columns under 
the torque requirements and under the momentum 
requirements apply to the condition wherein both 
the controlled and the residual modes are estimated. 
(Subsequent results in this paper pertain to this con- 
dition.) In the second columns, no estimates are used 
for any of these modes, and in the third columns, only 
the controlled modes are estimated. Naturally, the 
second columns for the flexible-mode requirements 


are blank because estimates of the controlled modes 
are required. 

When all of the. flexible modes are estimated, 
all cases shown in table X are stable. However, 
when certain modes are not estimated, some sensor 
arrangements produce unstable control. For rigid- 
body control, the lowest momentum requirements 
are obtained when none of the flexible modes are 
estimated, but an additional unstable case resulted. 
The case which gives the lowest momentum require- 
ment (M = 9712 ft-lb-sec) needs about 120 sec for 
the rigid-body mode ip to null. For all other cases 
shown in table X, the rigid-body modes null in 60 
sec. For flexible-mode control, the control require- 
ments are lowest when all the flexible modes are es- 
timated. Also, smoother control time histories are 
produced when all modes are estimated, as shown by 
the example in figure 5. 

In the bottom half of table X, the case giving 
the lowest requirement (M = 2602 ft-lb-sec) would 
seem to be the best for flexible-mode control, but 
figure 6 shows that this case results in poor modal 
response. Only the first and third modes are shown; 
however, a similar effect is produced in the second 
mode. The oscillations shown eventually die out 
because of the small amount of natural damping. 
In all other cases shown for flexible-mode control in 
table X, the flexible modes null within about 10 sec. 

For control of the rigid-body modes, the inclu- 
sion of flexible-mode initial disturbances produces 
only small increases in the momentum requirements 
(results in parentheses in table X). The maximum 
torque requirements, however, increase by about 
25 percent. 

Determination of Initial Estimates 

As previously stated, the results presented herein 
are determined with the assumption that the initial 
estimates of the modal displacements are 90 percent 
of the true values. Figures 7 and 8 show the length 
of time required for the observer (with controls off) 
to determine estimates to this accuracy. Figure 7 
shows time histories of the errors in the estimates for 
rigid-body control; that is, initial disturbances are 
included only in the rigid-body modes. Inspection 
of the figure shows that about 20 sec are required 
to obtain good estimates. Excessive errors can be 
noted for the residual-mode estimates; however, as 
previously stated, initial estimates for the residual 
modes are not necessary inasmuch as these modes 
are not in the control law. For flexible-mode control, 
good estimates are obtained after about 10 sec. (See 
fig- 8 .) 
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Example Time Histories for Rigid-Body 
Control 

Three-axis control. Time histories are given in 
figures 9 to 16 to compare the responses for the de- 
coupled and LQR procedures and to show the effects 
of inoperative control for these procedures. The re- 
sults pertain to rigid-body control and also illustrate 
the effects of lag in the system. As previously noted, 
the observer is in the loop for all cases. The initial 
conditions are given in the “General Considerations” 
section. Figure 9 contains the control responses for 
the two control procedures; lag is not included in the 
system. The modal responses resulting from the two 
control procedures are essentially identical; examples 
of these responses are shown in figure 10(c). Exam- 
ination of figure 9 shows that the LQR procedure 
basically relies on all six control actuators. For the 
decoupled procedure, the bulk of the control effort is 
performed by the three-axis CMG at the bottom of 
the column; hence, the CMG at the top may not be 
essential. 

Figure 10 shows the effect of assuming that the 
three- axis actuator at the top of the column is inoper- 
ative. This condition was simulated on the computer 
by zeroing out the first three rows in the feedback 
gain matrix. The other gains in the matrix were left 
unchanged; that is, the control law was not altered. 
As shown in figure 10, control of the system is still 
maintained by the one three-axis CMG, Comparison 
of the two cases for this CMG shows only small differ- 
ences in the torque time histories as illustrated by the 
small oscillations in two of the curves. The momen- 
tum time histories show that the one remaining CMG 
automatically compensates for the inoperative CMG, 
The rigid-body modal responses are essentially the 
same for both cases (fig. 10 (c)). The magnitudes 
of the flexible-mode transients for the inoperative- 
control case remain relatively low; however, the first 
flexible mode develops a slowly damped oscillation of 
very low magnitude. (As expected, when the three- 
axis actuator at the bottom of the column is assumed 
to fail, the system becomes unstable.) 

The results of incorporating a time lag r of 10 sec 
into the two rigid-body control procedures are pre- 
sented in figures 11 and 12. The lag is incorporated 
in the control input which enters both the plant and 
the estimator. Comparison of the time histories in 
figure 11 with those of the no-lag case (fig, 9) shows 
that the control requirements are considerably re- 
duced when lag is included. The rigid-body modal 
responses for the cases in figure 11 are essentially 
identical to those shown in figure 12 (c). 

As in figure 9, figure 11 indicates the CMG at 
the top of the antenna column is used sparingly for 


decoupled control. Hence, this case was investigated 
for the effect of an inoperative CMG at the top of the 
column. (See fig. 12 .) The simulation procedure was 
the same as that used for figure 10. The results in fig- 
ure 12 are similar to those in figure 10. However, note 
that the small oscillations in the torque time histories 
for the no-lag case (fig. 10 (a)) are eliminated by in- 
cluding lag. Comparing figure 12 (c) with figure 10 (c) 
shows that with lag included, more time is required 
to null the rigid-body responses. The comparison 
also shows that the magnitudes of the flexible-mode 
transients are substantially reduced when lag is in- 
cluded. The slowly damped oscillation in the first 
flexible mode (fig. 10 (c)) is virtually nonexistent in 
the case with lag. 

Single-axis control. Figures 13 to 17 are exam- 
ples of time histories relating to decoupled control 
of single-axis rigid-body modes. The control require- 
ments are summarized in table XI. In figure 13, the 
control responses are presented for control of the 
rigid-body mode 6, Results are similar for control 
of (j) and *0. The resulting 6 response for figure 13 is 
the same as that shown in figure 10(c). The other 
two rigid-body modes, as well as the flexible modes, 
are not materially affected. The data in figure 13 
correspond to the use of one three-axis attitude sen- 
sor at the top of the column. In the figure, the small 
torques for the z-sxis control are required because of 
the antenna product of inertia Ixz^ The normal con- 
trol arrangement for the two CMG’s was changed by 
moving the bottom CMG to location 4 . (See fig. 1 .) 
This arrangement produces much smaller effects in 
the second residual mode A 5 , as shown in figure 14. 

It should be noted that for single-axis rigid- 
body control, the residual modes were not estimated. 
Eliminating these estimates has little effect on the 
control requirements. (See table XL) Also, compar- 
ing figure 13(a) with figure 15 shows that the z-oxis 
control time histories are smoother when the residual 
modes are not estimated. (The effects are similar for 
the other ^-axis control actuator.) 

Control (about the x- and ^-axes) was maintained 
when the two z-exis control actuators for the case 
considered in figure 13 were assumed to be inopera- 
tive. This condition was simulated on the computer 
by zeroing out the third and sixth rows in the feed- 
back gain matrix. The other gains in the matrix were 
left unchanged; that is, the control law was not al- 
tered. (It should be noted that the system becomes 
unstable if only one z-axis actuator is inoperative.) 
The only noticeable effects of the inoperative actua- 
tors are illustrated in figure 16. As shown, a slowly 
damped oscillation of very small magnitude devel- 
oped for the first flexible mode Ax. In addition. 
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an extremely small steady-state offset results in the 
rigid-body mode 

The effect of eliminating the y- and z-sxis attitude 
sensors for the case considered in figure 13 was stud- 
ied. (See fig. 17.) Rigid-body control is maintained; 
however, the 0-axis control requirements increase, as 
shown by. the example in figure 17. The peak torques 
increase from about 4 to about 30 ft-lb, with a slowly 
damped oscillation. The overall control requirements 
increase by about 10 percent. The other effect no- 
ticed is in the first flexible mode, as shown by the 
oscillations in figure 17(b). As for the other flexible 
modes, A 2 is not affected and the small transient in 
A3 damps out rapidly. 

As previously mentioned, the summary in ta- 
ble XI shows that for single- axis rigid-body control, 
the control requirements are about the same whether 
or not the residual modes are estimated. Also, the 
control-actuator arrangement is not of importance, 
except for control of the rigid-body mode t/). For this 
case, the excessive requirements for locations 1 and 4 
were caused by large gains on 1/; in the feedback gain 
matrix F. (See table IV.) These gains resulted from 
the fact that the 0-mode slope values at locations 1 
and 4 for the first flexible mode A\ are about the 
same magnitude. 

Rigid-Body Command Requirements 

Decoupled-control requirements for 0.01-rad com- 
mands in the three rigid-body modes are presented 
in figure 18. The initial estimates are assumed to 
be 0.001 rad, or 10 percent of the commanded val- 
ues. For the computer simulations, it was necessary 
to scale the rigid-body command values v by multi- 
plying them by the squares of their respective closed- 
loop frequencies. 

Except where noted, the results in figure 18 corre- 
spond to the use of one three-axis attitude-rate sensor 
at location 1 and one three-axis attitude sensor at lo- 
cation 4. The single data point on figures 18(a) and 
18(b) for the one three-axis attitude sensor (at loca- 
tion 1), normally used for nulling rigid-body distur- 
bances, shows about the same control requirements; 
however, the control time histories are not as smooth. 
As noted in figures 18(a) and 18(b), various values of 
time lag were incorporated into the feedforward gain 
matrix G. (This matrix is shown in table V.) 

The curves in figures 18(a) and 18(b) depict the 
control requirements plotted against the amount of 
time lag used in the control feedback gains with 
TQ = 10 sec. When compared with the curves, the 
data points for both rp and tq = 0 sec illustrate 
the benefits of incorporating lag into the system. 
Also, the data points (open symbols) for tq = 0 
and tq ^ 5 sec indicate that feedforward lag is more 


effective than feedback lag. (For = 0 sec, large 
overshoot occurs in the rigid-body responses.) For 
reduced closed-loop frequencies, the solid symbols 
for TQ = 0 show a marked decrease in the control 
requirements and there is no overshoot; however, this 
condition requires instantaneous values for torque at 
^ = 0 sec. The solid symbols at = 6 sec show 
lower control requirements, but the time required to 
reach rigid-body steady state is doubled. 

The curves in figures 18(a) and 18(b) do not apply 
to Tp = 0. Data for this value are shown by single 
data points (for rp = 0, T(9 = 10) in these figures. 
The point in figure 18(b) shows a lower momentum 
requirement but, here again, this condition requires 
instantaneous values for torque at ^ = 0 sec. 

The curve in figure 18(a) shows that the torque 
|Tmax| is insensitive to rp over a wide range of values. 
This is not the case for momentum (fig. 18(b)), for 
which the best result occurs at about Tp = Q sec. 

Example Time Histories for Flexible-Mode 
Control 

Time histories are given in figures 19 to 23 to 
show various effects on the decoupled and LQR pro- 
cedures for flexible-mode control. The results per- 
tain to the effects of (1) initial estimate, (2) lag, and 
(3) inoperative controls. 

Effect of accuracy of initial estimate. Figures 19 
and 20 illustrate the effect of the initial estimates 
of the flexible-mode disturbances for the decoupled 
and LQR procedures, respectively. Parametric data 
on this effect are presented in a subsequent section 
which includes the effect on rigid-body control. In 
figures 19 and 20, comparison of the data for the 
two values of the initial estimate shows that perfect 
knowledge of the initial disturbances (initial estimate 
of 100 percent) results in the higher control require- 
ments for both control procedures. The differences 
in the torque and momentum requirements are not 
as great, however, for the LQR procedure as they 
are for the decoupled procedure. As shown in fig- 
ures 19(c) and 20(c), the flexible-mode responses are 
essentially unaflPected by the value of the initial es- 
timate. Noticeable differences are apparent in the 
rigid-body transient responses for decoupled control 
(fig. 19(c)). In figure 20(c), the rigid-body transients 
are relatively large for the LQR-control procedure, 
regardless of the value used for the initial estimate. 

Effect of lag. The effect of incorporating lag into 
the system for flexible-mode control can be seen by 
comparing figures 21 and 22. Comparisons of the 
two control procedures are also shown in each figure. 
Figure 21 is an example of a no-lag case, whereas for 
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figure 22 a time lag of 10 sec is included in the control 
feedback gains. 

Comparison of figures 21 and 22 shows that the 
lag reduces the maximum torque values by about 
50 percent. The momentum values are reduced by 
a lesser amount because of the additonal torque 
oscillations induced by the lag. 

Examination of figure 21(c) shows that the 
flexible-mode responses (no lag included) are about 
the same for the two control procedures. Although 
not shown, this result is the same for the case in 
which lag is included. In the lag case, the responses 
require a longer time to null than in the no-lag case. 
(An example of the responses is shown in fig. 23.) 
As shown in figure 22(a) for the actuator at the top 
of the antenna column, inclusion of lag resulted in 
virtually identical torque responses for the two con- 
trol procedures. Comparing the rigid-body responses 
in figure 21 (c) for the two control procedures shows 
relatively large transients for the LQR procedure. 

Effect of inoperative controls. Results in figures 21 
and 22 show that the decoupled procedure basically 
relies on all six control actuators, whereas, for the 
LQR procedure, the bulk of the control effort is 
performed by the three-axis CMC at the top of the 
antenna column. (The presence or absence of lag did 
not affect these results.) These results are opposite 
to those for rigid-body control with regard to both 
actuator location and control procedure. It is shown 
in figure 22 that for LQR control, lag virtually cancels 
the control effort for the three-axis CMC at the 
bottom of the antenna column. Figure 23 shows 
the effect of assuming that this CMG is inoperative 
when LQR control is used. The method used for 
this simulation was discussed in a previous section 
on rigid-body control. As shown in figure 23, control 
is still maintained by the one three-axis CMG, 
with only negligible changes in both the control 
requirements and modal responses. (As expected, 
when the three-axis actuator at the top of the column 
is assumed to fail, the system becomes unstable.) 

Ejffect of Control on Residual Modes 

Figure 24 shows the effect of decoupled control 
on the residual modes. The LQR-control results are 
similar in appearance to those shown for decoupled 
control, with one exception. For rigid-body control 
with no lag (compare fig. 24(b)), LQR control reduces 
the maximum peaks in A 4 and Ae by greater than 
50 percent and in A 5 by about 20 percent. 

Figure 24(a) pertains to flexible-mode control 
with no lag. Figures 24(b) and 24(c) pertain to 
rigid-body control with and without lag. Standard 


closed-loop dynamics were used for figure 24(a); LQR 
closed-loop dynamics were used for figures 24(b) and 
24(c). The results of figure 24 are essentially the 
same, regardless of the set of closed-loop dynamics 
used. 

Flexible-mode control. In figure 24(a), the initial 
disturbances are selected as 0.1, 0.5, and 0.1 in. for 
modes 4, 5, and 6 , respectively. As shown, these 
disturbances damp out because of the small values 
of natural damping in these modes (table I). The 
behavior of these modes is the same regardless of the 
magnitudes of the initial disturbances. Under certain 
conditions, such as type, number, and location of 
sensors and number of flexible modes estimated, one 
or more of the residual modes may become unstable, 
leading to an unstable control system. (See table X.) 
As far as flexible-mode control is concerned, the 
inclusion of lag in the system had little or no effect 
on the residual-mode time histories. 

Rigid-body control. The results in figures 24(b) 
and 24(c) apply to either rigid-body commands or 
nulling rigid-body disturbances. No initial distur- 
bances are assumed to exist in any of the flexible 
modes, including the residual modes. (Similar re- 
sults are obtained even when initial disturbances are 
present.) Although the disturbances damp out for 
the no-lag case, it is shown in figure 24(b) that rel- 
atively large amplitudes occur initially in the sec- 
ond residual mode (mode 5). With lag included 
(fig. 24(c)), the maximum amplitudes in all three 
residual modes are reduced by an order of magni- 
tude. As mentioned above, under certain conditions 
one or more of the residual modes may become un- 
stable. (See table X.) 

Parametric Summary of Control 
Requirements 

Effect of initial estimate. Some time-history 
examples were previously presented (figs. 19 and 20 ) 
on the effect of the accuracy of the initial estimates. 
Figure 25 presents results of control requirements 
plotted against initial estimates. The rigid^body 
results (fig. 25(a)) are shown for decoupled control 
only; results are similar for LQR control. For rigid- 
body control, the lowest torque requirements occur 
at an initial-estimate value of about 90 percent. Over 
the range of initial estimates, the iTmaxI values are 
produced, for the most part, by the a:-axis control 
actuator (at location 2 ) at a time of about 2 sec. 
(For initial estimates of 95 percent and above, the y- 
axis actuator produces the maximum torque.) As for 
momentum, figure 25(a) shows the best results for an 
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overestimate of about 10 percent. (This result would 
change somewhat by altering the estimator feedback 
gains.) The overestimate could be determined by 
taking the value of a “good” estimate, say 85 to 
90 percent, and adding 20 percent. A good estimate 
could be based on a priori information such as that 
shown in figures 7 and 8. A good estimate could 
also be established by monitoring the estimator error 
£ = y — C'ic. 

The data for fiexible-mode control (fig. 25(b)) 
are shown for the three conditions noted in the 
key. All data for |Tmaxl up to an initial estimate 
of about 100 percent apply to the «-axis control 
actuator at location 1. The iTmaxI values are about 
10 percent higher for decoupled control with the 
standard closed-loop dynamics. Also, iTmaxI for 
decoupled control is insensitive to the initial estimate 
over the range up to about 100 percent. The curve 
for LQR control does not vary appreciably, with the 
lowest value of |Tmax| occurring at an initial estimate 
of 75 percent. 

The lowest momentum requirements occur at ini- 
tial estimates of 50 to 70 percent, depending upon 
the control procedure used. Also, for initial estimates 
from 0 to 100 percent, the LQR-control momentum 
requirements are 3 to 15 percent lower than those for 
decoupled control. 

Effect of lag. Figure 26 presents results show- 
ing the effect of lag on control requirements. (As 
previously noted, the lag was incorporated into the 
control input which enters both the plant and the es- 
timator.) The data are shown for decoupled control 
only in figure 26(a); LQR-control momentum results 
are about the same but |Tmaxl results are lower than 
decoupled-control results. (See table XII.) As a mat- 
ter of interest, results are also included for the ideal 
case (no estimator). Perfect knowledge of the state 
variables for feedback was assumed and no residual 
modes were included in the model. 

Results for fiexible-mode control are presented for 
the two control procedures in figure 26(b). Standard 
closed-loop dynamics were used for the decoupled 
control; results with LQR closed-loop values were 
about the same. (See case 2C in table XII.) Results 
for the ideal case were the same as those shown, 
except that the decoupled-control requirements were 
several percent higher at lag values of 0 < t < 2 sec. 

The data in figure 26(b) pertain to use of initial 
estimates of 90 percent. The value used for the initial 
estimate is not crucial, however, when lag is included 
in the system. For example, with a lag of 10 sec, a 
value of 60 percent for the initial estimate results in 
approximately the same magnitudes as those shown 
for both torque and momentum. 


In comparing the two control procedures in fig- 
ure 26(b), note the pronounced decrease in momen- 
tum for the LQR-control procedure as t varies from 
0 to 5 sec. This effect of lag is similar to the effect of 
decreasing the flexible-mode closed-loop damping ra- 
tios for this procedure, as previously discussed for ta- 
ble VIII. The disadvantage of the LQR-control proce- 
dure, as previously discussed regarding the transient 
responses, is evident when lag is included. Table XIII 
(case 2C) shows rigid-body transients of about one 
order of magnitude higher than those for decou- 
pled control. Table XIII also shows that for rigid- 
body control, the flexible-mode transients are one to 
two orders of magnitude higher for the LQR control 
procedure. 

In table XII, it is shown that for rigid-body con- 
trol, use of the two three-axis sensors produces the 
lowest control requirements. Case 2D is included 
in table XII to indicate the lower control require- 
ments for LQR flexible-mode control when only the 
one three-axis CMC (at location 1) is in operation. 
(See also fig. 23.) Table XIII, however, shows that 
for this case the rigid-body transients are close to 
one order of magnitude higher. (Compare cases 2D 
and 2C.) 

Concluding Remarks 

An analysis has been performed to compare de- 
coupled and linear quadratic regulator (LQR) proce- 
dures for the control of a large, flexible space an- 
tenna. A state estimator is used in the loop for 
both procedures. Control involved (1) commanding 
changes in the rigid-body modes, (2) nulling initial 
disturbances in the rigid-body modes, or (3) nulling 
initial disturbances in the first three flexible modes. 
Control was achieved with two three-axis control- 
moment gyros (CMG’s) located on the antenna col- 
umn. The results illustrated various effects on con- 
trol requirements for the two procedures. These 
effects included errors in the initial estimates of state 
variable, variations in the type, number, and loca- 
tion of sensors, and deletions of state-variable esti- 
mates for certain flexible modes after control activa- 
tion. The advantages of incorporating a time lag in 
the control feedback were also illustrated. In addi- 
tion, the effects of inoperative-control situations were 
analyzed with regard to control requirements and re- 
sultant modal responses. Comparisons were included 
which showed the effects of perfect state feedback 
with no residual modes (ideal case). Time-history 
responses were presented to illustrate the various ef- 
fects on the control procedures. Results of the study 
indicate the following: 

1. Compared with the ideal case (perfect esti- 
mates of state variables and no residual modes), 
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control requirements were appreciably higher when 
an observer was incorporated in rigid-body control. 
Further, for flexible-mode control, the requirements 
were generally somewhat lower than for rigid-body 
control. 

2. Accuracy of the initial estimates of the state 
variables had a significant effect on control require- 
ments, especially for rigid-body control. In fact, an 
overestimate of about 10 percent for rigid-body con- 
trol produced the lowest momentum requirements. 
For flexible-mode control, initial estimates of from 
50 to 75 percent resulted in the lowest control re- 
quirements. About 10 sec were required before con- 
trol activation to obtain good estimates for flexible- 
mode control, and 20 sec were required for rigid-body 
control. 

3. For flexible-mode control, the LQR-control re- 
quirements were less than those for decoupled con- 
trol, with a pronounced reduction for low values of 
closed-loop damping. For rigid-body control, the 
momentum requirements were about the same for 
both control procedures; however, the LQR maxi- 
mum torque requirements were much less. 

4. Linear quadratic regulator (LQR) control had 
more of an effect on the transients which occur in the 
initially undisturbed modes. For either rigid-body 
control or flexible-mode control, the LQR transients 
were generally an order of magnitude larger than 
those for decoupled control. Incorporation of lag 
into the system can be used to reduce the transients 
considerably. 

5. The type, number, and location of sensors had 
only a small effect on control requirements when the 
control system was stable. Rate sensors cannot be 
used alone because the system is unobservable with- 
out attitude sensors. Two three-axis sensors (one for 
attitude and one for attitude rate) produced lower 
control requirements than one three-axis attitude 
sensor. For single-axis rigid-body control, a one-axis 
attitude sensor was adequate; however, a three-axis 
attitude sensor produced better results. (In practice, 
three-axis control is usually required.) 

6. With certain sensor arrangements, instability 
occurred when the residual modes were not esti- 
mated after control activation. This was also true 


when none of the flexible modes were estimated. For 
flexible-mode control, the torque and momentum re- 
quirements were generally lower when all the flexible 
modes were estimated. 

7. Under certain conditions, for either control 
procedure, one three-axis CMG could be made in- 
operative without materially changing the control re- 
quirements or the resultant modal responses from the 
normal case. 

8. A time lag incorporated into the system led 
to large reductions in the control requirements. The 
reduction was pronounced in the LQR momentum re- 
quirements for flexible-mode control. Incorporating 
lag in the control input automatically included lag 
in the estimator. For rigid-body control, lag was also 
effective in reducing the residual-mode amplitudes by 
one order of magnitude. 

9. For either rigid-body or flexible-mode control 
of initial disturbances, the control requirements de- 
creased as lag increased. For rigid-body commands, 
the requirements were essentially invariant with lag 
values of up to about 10 sec. For command control, 
it was essential that lag be included in both the feed- 
forward and feedback control gains. 


NASA Langley Research Center 

Hampton, VA 23665-5225 

June 13, 1985 

References 

1. Sullivan, Marvin R.: LSST (Hoop/Column) Maypole 
Antenna Development Program. NASA CR-3558, Part 2, 
1982. 

2. Young, John W.; Hamer, Harold A.; and Johnson, 
Katherine G.: Decoupled- Control Analysis of a Large 
Flexible Space Antenna With Linear Quadratic Regulator 
Comparisons. NASA TP-2293, 1984. 

3. Joshi, Suresh M.: Control Systems Synthesis for a 

Large Flexible Space Antenna. Acta Astronaut, vol. 10, 
no. 5-6, May- June 1983, pp. 365-380. 

4. Kwakernaak, Huibert; and Sivan, Raphael: Linear Op- 
timal Control Systems. John Wiley & Sons, Inc., c.1972. 

5. Armstrong, Ernest S.: ORACLS — A System for Linear- 
Quadratic- Gaussian Control Law Design. NASA TP-1106, 
1978. 


13 



TABLE I. ANTENNA CHARACTERISTICS 

Ix = 1.9577 X 10^“ IbV; ly = 1.9657 x 10^® Ib-in^;' 
h = 1.499 X 10^^^ Ib-in^; = 0.8357 x 10® Ib-in^; 
w = 10020.3 lb 


Mode 

ojn, rad/sec 

T^, SGC 


run, Ib-sec^/in. 

Description 

1 

0.7466 

8.42 

0.01 

153.157 

First torsion 

2 

1.3460 

4.67 

.01 

5.233 

First bending, 
x-z plane 

3 

1.7025 

3.69 

.01 

3.073 

First bending, 
y-z plane 

4 

3.1813 

1.98 

.01 

.305 

Surface torsion 

5 

4.5294 

1.39 

.02 

1.993 

Second column bending, 
y-z plane 

6 

5.5905 

1.12 

.02 

723.522 

Surface-hoop torsion, 
feed mast torsion 
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TABLE II. B MATRICES USED IN CURRENT ANALYSIS 


(a) Three-axis CMG’s at locations 1 and 2 

Bi8x6 


0. 

9. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0. 

•1^737917f-07 

0. 

.110039876-09 

.197379176-07 

0. 

.110039876-09 

0, 

.19658120E-07 

0 . 

0. 

.196581206-07 

0. 

.110039S7E-OQ 

0. 

.257777986-07 

.110039876-09 

0. 

.257777986-07 

.l^712335E-0? 

0, 

-.117454^16-03 

-.423205426-06 

0. 

-.377761956-06 

C, 

-.3^6810036-03 

0. 

0 . 

.318386146-04 

0. 

•61143590F-03 

0. 

.897564476-04 

-.618855136-04 

0. 

.200983116-05 

-.3c»176796F~06 

0. 

.441990446-05 

.887460346-07 

0. 

-.102762376-03 

-.l?69352ir-03 

0. 

-•114521516-04 

-.219815676-02 

0. 

-.835529376-06 

• 5958a099C-07 

0. 

-.342740296-05 

.716896366-07 

0. 

-.135514966-03 




(b) Three-axis CMG’s 

at locations 1 and 4 





Bi8x6 



0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

.197379176-07 

0. 

•110039876-09 

•197379176-07 

0. 

.110039876-09 

0. 

.196581206-07 

0. 

0. 

.196501206-07 

0. 

.110039876-09 

0. 

.257777986-07 

.110039876-09 

0. 

.257777986-07 

.147123256-05 

0. 

-.117454416-03 

-.771630216-06 

0. 

-.715602196-04 

0 • 

-.366B10C3E-03 

0. 

0. 

.103343246-03 

0. 

.611435006-03 

0. 

.897564476-04 

-.164736356-03 

0. 

.555204306-04 

-.391767986-06 

0. 

.441990446-05 

.785225476-07 

0. 

.121484186-03 

-.126085216-03 

0. 

-.114521516-04 

.154323316-03 

0. 

-.742944766-05 

.595890996-07 

0. 

-.342740296-05 

-.594732216-bB 

0. 

-.339336246-04 








TABLE IV. F MATRICES FOR DECOUPLED CONTROL 
[Standard closed-loop dynamics used] 


(a) Three-axis CMG’s at locations 1 and 2 

^6x18 


j-.<.671?'V?7f + C5 
Row no. 1\0* 

(-'.96l93614r + 03 

jo. 

2 < 0 . 

io. 

( .10685980F+04 
3<0. 

( .73817760E+04 
(-.45993871E+06 

V .9619361^E+03 

! 0. 

C. 

0 . 

.10941967F+04 

0. 

-.73817760E+04 


0 . 

0 . 

0, 

-.40627667E+05 

0 . 

.22234285E+04 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

-♦46306799E+06 

0 . 

-.22264285E+04 

0 . 

0 . 

0 . 


.11921449E+04 

0 . 

-.11606671E+04 

0 . 

0 . 

0 . 

.12632075E+04 

0. 

-.18780960E+02 

.97063982E+03 

0 . 

.11606671E+04 

0 . 

0 . 

0 . 

-•38920317E+C6 

0 . 

.18780960E+02 


-.84082369E+06 

0 . 

0 . 

0. 

0 . 

0. 

.19234584E+05 

0 . 

0 . 

-.827689638+07 

0 . 

0 . 

0 . 

0. 

0 . 

.19695540E+05 

0 . 


0 . 

0 . 

0. 

0 . 

-.73129801E+06 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.84252233E+07 

0 . 

0 . 

0 . 

0 . 


0 . 

.21458608E+05 

0. 

0 . 

0 . 

0. 

0. 

.22737735E+05 

0. 

0 . 

.17471517E+05 

0 . 

0 . 

0 . 

0. 

0 . 

-.70056570E+07 

0 . 


(b) Three-axis CMG’s at locations 1 and 4 

^6x18 


Row no 


-.107682?3E+0f 
7383E+03 


(-.1076 
. 1 < 0 . 

(-.830*^ 

1^: 

}“• 

I”: 

li; 

( .22701 

Jo. 

V-.1R«30 


-.11632138E+03 

<^«30121F+0*^ 
39896 «^91f+06 

830*>7363P + 03 


227010616+04 
30121F+05 


0 . 

0 . 

0 . 

-.11181I514E + 06 

0* 

.1889^1158+04 

0 . 

0. 

0. 

0. 

0. 

0. 

-.3968P051E+06 

0. 

-.1B895115F+04 

0. 

0. 

0. 


.15260150E+04 

0. 

-•10068196E+04 

0. 

0. 

0. 

.604O3111E+06 

0. 

-.49203577E+02 

.636769676+03 

0. 

.10068196E+04 

0. 

0. 

0. 

-. 992871078+06 

0. 

.492035778+02 


-.193828016+07 

0. 

0. 

0* 

0. 

0 . 

-.209378496+04 

0. 

0. 

-.718144046+07 

0. 

0. 

0. 

0 . 

0* 

,4X0239006+05 

0, 


0 . 

0. 

0 . 

0. 

-.201267268+07 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.714384026+07 

0 . 

0 . 

0 . 

0 . 


.274682708+05 

0 . 

0. 

0 . 

0 . 

0. 

.108887608+08 

0 . 

0 . 

.114618548+05 

0. 

0 . 

0 . 

0 . 

0 . 

-.178716798+08 

0 . 



00 


•46712427F+0T 

0 . 

-.10685BRQE+06 

.45993B71E+0P 

C. 

-.10941967E+06 


TABLE V. G MATRIX USED FOR RIGID-BODY COMMANDS 
[Three-axis CMG’s at locations 1 and 2] 

Gex6 


0. 

.40627667E+07 

0 . 

0 . 


.463O6799E+0e 


0 . 


-.119214496+06 

0 . 

-.12632075E+06 

-.97063982E+05 

0 . 

.38920317E+08 


•11107667F+0A 

0 . 

-•8f234396E+04 

~.lllC7667E+04 

0 . 

• 85234396E^-04 


-.25084746E+04 

0 . 

0 . 

.25084746E+04 


0 . 


•14B20496E+04 

0 . 

.23981307F+02 

-.148204968+04 

0 . 

-.23981307E+02 



TABLE VI. F AND K MATRICES FOR DECOUPLED- AND LQR-CONTROL COMPARISONS FOR CASE 2 OF TABLE III 

[LQR closed-loop dynamics used] 

Fexis 


Row no. 



4 + 

PI + 


0367'^('PP+OA 

6^51.92^1F404 


1474S42F40^ 


0, 

0, 

0 . 

'I-; 


1 Of'l 70F + 0^ 

A7«:i Q'JPI F+04 


0 . 

0 . 

0 . 

f O'iQOOT^F + O'^ 

0 . 

,3B«27P20F+C4 

0 . 

0 . 

0 . 

0 . 

C. 

0 . 

BF + 06 

0 . 

?3‘^27fl?o£ + G4 


0. 

0. 

0 . 


.1 4Q0090BF+04 

.21542Q86F+0T 

0. 

.136144196+02 

0. 

-.B31524ieP+06 

0. 

.240200456+05 

-.24803124P+04 

0. 

0. 

0. 

0. 

0. 

-.809976426+01 

0. 

0. 

0. 

-.8164Q130F+06 

0. 

0. 

0. 

0. 

0. 

.15789134F+C4 

-.16530^*?5F + 0? 

0. 

.220297326+00 

0 • 

.190?lfi4oc + 0‘^ 

0. 

.254518586+05 

-.401343731^ + 02 

0. 

0. 

0. 

21322606 + 04 

-.21542984P+01 

0. 

-.136144196+02 

0. 

-.818737;>2r + 07 

0. 

.195570306+05 

.248031246+04 

0. 

0. 

0. 

0. 

0. 

. R0997642E+01 

0. 

0. 

0. 

-.940672R6F+C7 

0. 

0. 

0. 

0. 

0. 

-.486474336+06 

.16530”55F+07 

0. 

-.220297326+00 

0. 

.lQ4777f>PF+05 

0. 

-.784189746+07 

.401343736+02 

0. 

0. 

0. 


Kexis 


(O 


Row no. 


1 

2 

3 

4 

5 

6 



29695400F406 

37985200F+03 

64377600F+04 

61915100E+04 

33254100E+06 

80414800E+02 

12852600E+04 

935697006+02 


0. 

0 . 

0 . 

-.293977006+06 

0 . 

.363859006+04 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.337012006+06 

0 . 

-.342500006+03 

0 . 

0 . 

0 . 


.549278006+04 

0 . 

-.265962006+04 


0 . 

0 . 

0 . 


-.164308006+03 

-.496916006+07 

0 . 

0 . 

0 . 

0 . 


0 . 

0 . 

0 . 


.298704006+03 

-.493832006+07 


0 . 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


166435006+06 

324594006+03 

921316006+03 

281736006+03 

205592006+06 

362612006+01 


.446188006+03 

-.111494006+06 

0 . 

.915410006+00 

-.535499006+07 

0 . 

0. 

0 . 

0 . 

.519965006+03 

-.124421006+05 

0 . 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

.319770006+03 

-.543606006+07 

0 . 

0 . 

0 . 

0 . 


-.217185006+03 

.105906006+06 

0 . 

0 . 

0 . 

0 . 


-.160382006+03 

-.327119006+07 

0 . 

-.234356006+03 

-.134411006+05 

0 . 

0 . 

0 . 

0 . 

-.221741006+02 

-.380137006+07 

0 . 
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TABLE VII. C MATRIX FOR THREE-AXIS ATTITUDE-RATE SENSOR AT LOCATION 1 
AND THREE-AXIS ATTITUDE SENSOR AT LOCATION 4 

C6xl8 


10 . 

Row no . 1 \ 0 • 

{ ♦22532850f-03 


. 17983610F-01 

• lOOOOOOOE + 01 

• 23921 ^70E-07 


6< •37009500E-04 

(o. 


0 . 

0 . 

0 . 

0 . 

0* 

-.19195<rl0E-02 

0* 

0 . 

0 . 

0 . 

.30756450E-03 

0, 

.lOOOOOOOE+01 

0 . 

0, 

0 . 

-.14806800E-04 

0 , 


♦ l«'790460E-02 

0 . 

0 . 

0, 

0* 

0* 

*27583190F-03 

0 . 

-•43030160E-05 

0, 

0 . 

0 . 

0 . 

.lOOOOOOOE+01 

-.245'^1710E-01 

0 . 


0 . 

.lOOOOOOOE+01 

-♦11935500E-06 

0* 

0* 

0 . 

0* 

0, 

.13464720E-05 

-.11816080E-03 

0 . 


►10959970E-01 


0, 

0 . 

-.25308020E-03 

0* 

. lOOOOOOOE+01 

0 . 

0, 

0, 

-.22823800F-04 

0 . 

0, 

0* 

.54079040E-03 

0, 

0* 

0. 

0. 

0. 


.43113600£*04 


•lOOOOOOOE+01 

-.24797850E-02 

►•50625030E-03 


.17061950E-03 



TABLE VIII. DECOUPLED- AND LQR-CONTROL REQUIREMENTS WITH AND WITHOUT OBSERVER 
[Numbers in parentheses indicate control with observer; residuals estimated] 


Case 

No. of 
three-axis 
sensors 

Initial disturbances of — 

LQR-control 
requirements for — 

Decoupled-control 
requirements for — 

Rigid-body 
modes, rad 

Flexible modes, 
in. 

|Trnax|) 

ft-lb 

Momentum, 

ft-lb-sec 

|7max|) 

ft-lb 

Momentum, 

ft-lb-sec 

1 

1 

0.01 

0 

292 

9895 

488 

9910 





(318) 

(11 736) 

(509) 

(12056) 


2 

0 

1 

276 

3543 

278 

3850 





(258) 

(3440) 

(270) 

(3 678) 

2 

1 

0.01 

0 

281 

9945 

486 

9876 





(301) 

(11 768) 

(501) 

(11 798) 


2 

n 0 

1 

210 

2400 

212 

2850 





(202) 

(2373) 

(209) 

(2 750) 


1 

0 

1 

(195) 

(2228) 

(230) 

r (3091) 

3 

1 

0.01 

0 

267 

9976 

485 

9860 





(276) 

(11707) 

(495) 

(11555) 


2 

0 

1 

63 

1355 

65 

2 378 





(68) 

(1367) 

(61) 

(2312) 




TABLE DC. MAXIMUM TRANSIENTS DURING DECOUPLED AND LQR CONTROL WITH OBSERVER IN SYSTEM 


Case 

No. of 
three-axis 
sensors 

Initial disturbances of — 

Maximum transients 
for LQR control, 
rad or in. 
(absolute) 

Maximum transients 
for decoupled control, 
rad or in. 
(absolute) 

Rigid-body 
modes, rad 

Flexible modes, 
in. 

1 

1 

0.01 

0 

Ai = 0.16 

Ai = 0.047 





A2 = 0.37 

A2 = 0.164 





A3 = 0.50 

A3 = 0.170 


2 

0 

1 

0 = 0.27 X 10“^ 

0 = 0.10 X lO"^ 





(f) = 0.35 X 10“3 

<t> = 0.35 X 10“^ 





Ip = 0.91 X 10~® 

i) = 0.42 X 10"^ 

2 

1 

0.01 

0 

Ai = 0.25 

Ai = 0.036 





A2 = 0.54 

A2 = 0.126 





A3 = 0.60 

A3 = 0.127 


2 

0 

1 

0 = 0.16 X 10“^ 

9 = 0.10 X 10-^ 





(p = 0.23 X 10“2 

<p = 0.35 X 10-^ 





= 0.58 X 10“2 

Ip = 0.42 X 10“^ 


1 

0 

1 

e = 0.15 X 10“^ 

9 = 0.40 X lO"'’ 





<p = 0.19 X 10“3 

(p = 0.70 X 10-3 





= 0.50 X 10“2 

= 0.76 X 10-^ 

3 

1 

0.01 

0 

Ai = 0.47 

Ax = 0.011 





A2 = 0.94 

A2 = 0.058 





A 3 = 1.06 

A 3 = 0.055 


2 

0 

1 

9 = 0.40 X 10“^ 

0 = 0.10 X 10-^ 





<p = 0.76 X 10“^ 

<p = 0.35 X 10-^ 





xp = 0.15 X 10-3 

Ip = 0.42 X 10-^ 




TABLE X. DECOUPLED-CONTROL REQUIREMENTS FOR VARIOUS COMBINATIONS OF SENSORS 
WITH AND WITHOUT ESTIMATES IN FLEXIBLE MODES 

[U indicates unstable; numbers in parentheses indicate initial disturbances included in flexible modes] 



Three-axis sensors 

|Tmax|, ft-lb, determined — 

Momentum, ft-lb-sec, determined — 




With 

Without estimates 

Without estimates 

With 

Without estimates 

Without estimates 

Initial 



estimates in 

in flexible 

in residual 

estimates in 

in flexible 

in residual 

disturbances 

Type 

Location 

flexible modes 

modes 

modes 

flexible modes 

modes 

modes 

= 0.01 rad 

Attitude 

1 

418 

443 

418 

11033 

10546 

10984 

(An = 0) 

Attitude 

2 

436 

U 

U 

11094 

U 

U 


Attitude 

4 

429 

U 

447 

11103 

U 

10890 


Attitude 

1, 2 

511 

u 

U 

10773 

u 

U 


Attitude 

1, 4 

475 

480 

475 

10762 

10407 

10845 


Attitude 

2,4 

470 

U 

U 

10765 

U 

U 


Attitude, 
attitude rate 

2, 1 

418 

u 

U 

10456 

u 

u 


Attitude, 
attitude rate 

4, 1 

417 (529) 

414 

429 (540) 

10438 (10505) 

9712 

9916 (10667) 

Ai, ^ 2,^3 = 1.0 in.; 

Attitude 

1 

247 


246 

3030 


3066 

A 4 , Ae = 0.1 in.; 

Attitude 

2 

219 


U 

2602 


U 

A 5 = 0.5 in. 

Attitude 

4 

241 


275 

2890 


3570 

^ = 0) 

Attitude 

1, 2 

237 


U 

2844 


U 


Attitude 

1, 4 

241 


246 

2917 


2977 


Attitude 

2,4 

228 

! 


U 

2 767 


U 


Attitude, 
attitude rate 

2, 1 

226 


U 

2 700 


U 


Attitude, 
attitude rate 

4, 1 

227 


243 

2716 


2919 


09 




TABLE XI. CONTROL REQUIREMENTS FOR CONTROLLING SINGLE-AXIS RIGID-BODY MODES 


Initial 

disturbances 

Control 

locations 

Residuals 

estimated 

l^maxl, 

ft-lb 

Momentum, 

ft-lb-sec 

9 = 0.01 

1,4 

No 

362 

3905 


1,4 

Yes 

363 

3907 


1, 2 

No 

395 

3848 


1, 2 

Yes 

418 

3868 

1— 1 

O 

o 

II 

1,4 

No 

354 

3880 


1,4 

Yes 

354 

3 880 


1,2 

No 

418 

3888 


1,2 

Yes 

418 

3888 

II 

o 

o 

1— ^ 

1,4 

No 

767 

13158 


1,4 

Yes 

765 

13209 


1,2 

No 

300 

3343 


1,2 

Yes 

300 

3371 


24 




TABLE XII. DECOUPLED- AND LQR-CONTROL REQUIREMENTS FOR CASE 2 WITH OBSERVER AND LAG 

[r = 10 sec; LQR closed-loop dynamics used] 


Case 

No. of 
three-axis 
sensors 

Initial disturbances of — 

LQR-control 
requirements for — 

Decoupled-control 
requirements for — 

Rigid-body 
modes, rad 

Flexible modes, 
in. 

l^maxh 

ft-lb 

Momentum, 

ft-lb-sec 

|7max|, 

ft-lb 

Momentum, 

ft-lb-sec 

2A 

1 

0.01 

0 

108 

9070 

184 

9164 

2B 

2 ^ 

0.01 

0 

96 

8894 

169 

8948 

2C 

2 

0 

1 

91 

1401 

96 

2470 

2D 

2* 

0 

1 

90 

1292 




*One three-axis control actuator inoperative. 


to 



TABLE XIII. MAXIMUM TRANSIENTS DURING DECOUPLED AND LQR CONTROL 
FOR CASE 2 WITH OBSERVER AND LAG 

[r = 10 sec; LQR closed-loop dynamics used] 


Case 

No. of 
three-axis 
sensors 

Initial disturbances of — 

Maximum transients 
for LQR control, 
rad or in. 
(absolute) 

Maximum transients 
for decoupled control, 
rad or in. 
(absolute) 

Rigid-body modes, 
rad 

Flexible modes, 
in. 

2A 

1 

0.01 

0 

Ai = 0.16 
A2 = 0.27 
As = 0.28 

Ai = 0.018 

A2 = 0.075 
As = 0.060 

2B 

2 

0.01 

0 

Ai = 0.19 
A2 = 0.22 

As = 0.22 

Ai = 0.004 
A2 = 0.010 
As = 0.007 

2C 

2 

0 

1 

0 = 0.17 X 10“^ 
(t> = 0.24 X 10-4 
^ = 0.16 X 10-3 

0 ^ 0.22 X 10-3 
0 = 0.79 X 10-3 
rp = 0.17 X 10-4 

2D 

2* 

0 

1 

0 = 0.90 X 10-4 
<t> = 0.15 X 10-3 
ij) = 0.50 X 10-3 



*One three-axis control actuator inoperative. 










Time, sec 

(b) With lag; r == 10 sec. 

Figure 3. Flexible-mode transient response for decoupled and LQR control of rigid-body modes with observer 
in system. 
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Time, sec 


(b) With lag; r = 10 sec. 


Figure 4. Rigid-body-mode transient response for decoupled and LQR control of flexible modes with observer 
in system. 

















ft-lb 



Time, sec Time, sec 

(a) Torque response. 

Figure 9. Decoupled- and LQR-control responses for rigid-body control with no lag. (LQR closed-loop dynamics 






ft-lb-sec Myi, ft-lb-sec M,,i, ft-lb-sec 


Two 3-axis actuators 


One 3-axis actuator 




Time, sec Time, sec 


(b) Momentum required for decoupled control. 
Figure 10. Continued. 






Decoupled 





Time, sec 


20 40 60 

Time, sec 


(a) Torque response. 

Figure 11. Decoupled- and LQR-control responses for rigid-body control with lag. (LQR closed-loop dynamics 




Two 3-axis actuators 


One 3-axis actuator 



(a) Torque response for decoupled control. 

Figure 12. Effect of inoperative CMG at top of antenna column for rigid-body control with lag. (LQR closed- 
loop dynamics used.) 
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ft-lb-sec My i, ft-lb“sec ft-lb-sec 


Two 3-axis actuators 


One 3-axis actuator 







(b) Momentum required for decoupled control. 
Figure 12. Continued. 
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in 


l6 


Time, sec 

(a) z-axis torque response at location 2. 



Time, sec 


(b) First-flexible-mode response. 

Figure 17. Effects of inoperative y- and z-axis attitude sensors for single-axis rigid-body decoupled control. 


(a) I 


Figure 18. Decoupled-control requiremeni 
represent 50-percent reduction in close 










Initial estimate = 70^ 

Initial estimate = 100^ 



Time, sec Time, sec 

(a) Torque response. 

Figure 19. Effect of accuracy of initial estimate for flexible-mode decoupled control. (Standard closed-loop 
dynamics used.) 
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ate = -jot 
ate = 100% 
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0 

800 
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/[omentum required. 




19. Continued. 





















Time, sec 


Time, sec 


(a) Torque response. 

Figure 21. Decoupled and LQR control responses for flexible-mode control with no lag. (LQR closed-loop 
dynamics used.) 
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LQR 

8001 
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jT 

<1> 
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A 
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w 600 


w 600 














-sec 


Decoupled 

■LQR 


600 

O 

<u 

W 450 


w 450 


w 450 


Time, sec 

Time, sec 

(b) Momentum required. 
Figure 22. Concluded. 







One 3-axis actuator 


Two 3-axls actuators 



Time, sec Time, sec 

(a) Torque response for LQR control. 

Figure 23. Effect of inoperative CMG at bottom of antenna column for flexible-mode control. (Conditions 
correspond to those of fig. 22.) 
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Time, sec 


Time, sec 


Time, sec 


(a) Flexible-mode control; r = 0. 


(b) Rigid-body control; r = 0. 


(c) Rigid-body control; r = 


Figure 24. Effect of decoupled control on residual modes. 
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